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ABSTRACT 

We study the sawtooth lattice of a coupled spin 1/2 Heisenberg system, a variant of 
the railroad trestle lattice. The ground state of this system is two-fold degenerate with 
periodic boundary conditions and supports kink antikink excitations, which are distinct in 
this case, unlike the railroad trestle lattice. The resulting low temperature thermodynamics 
is compared with the recently discovered Delafossites YC11O2.5 ■ 



1 



I. INTRODUCTION 



Several novel quantum toplogical excitations have been explored in the quantum an- 
tiferromagnetic spin systems in the last decade, expanding on the Anderson-Kubo spin 
waves. These include "spinons", i.e. spin 1/2 objects of Faddeev and Takhtajan [1], a name 
coined by Anderson, describing the excitations of the isotropic one-dimensional Heisenberg 
spin 1/2 antiferromagnet (AFM), with a concomitant two-parameter elementary excitation 
spectrum with an asymptotic four-fold degeneracy. The anisotropic Heisenberg AFM in 
the Neel ordered, massive phase contains "domain walls" found by Johnson, Krinsky and 
McCoy [2] from Baxter's solution of the XYZ model. These domain walls separate two 
regions of Ising like ordered states, and propagate as dressed fermions [3] (again these ex- 
hibit a four- fold degeneracy). The domain walls would broaden out and indeed the width 
would diverge with the correlation length as the isotropic point is reached, so the limiting 
excitation would have a delocalized character. It may be tempting to view the "spinons" 
as limiting cases of the domain walls, although the AFM order vanishes at the isotropic 
point. Yet another class of excitations were introduced by Shastry and Sutherland [SS] 
[4] in the context of models with broken translation symmetry [5]. These are topological 
quantum solitons separating two regions of broken translational symmetry and have again 
a four-fold degeneracy of two-parameter excitations. In addition, these are fairly compact 
objects with a width consisting of a few lattice constants. In view of the theoretical interest 
in these constructs, it is interesting to ask for experimental realizations of such systems. 

Recently a new family of Delafossite compounds have been synthesized which seem to 
be promising from this viewpoint. The YCuC>2 lattice consists of planes of coupled YO2 
octahedra linked by two-fold-coordinated bridging Cu + ions, which form a triangular 
planar array. It is possible to intercalate 2 ~ ions into the Cu + planes, forming different 
lattice symmetries depending on the amount of intercalant. For compositions up to and 
including YC11O2.5 > the planar 2 ~ form an orthorhombic structure [6]. At the upper 
limit of this range one has a magnetic insulator, with all of the copper ions converted to 
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Cu 2+ . For compositions with additional oxygen beyond O2.5 , the structural symmetry 
becomes hexagonal and mobile hole carriers appear, rendering the CuO x planes weakly 
conducting [7]. 

Preliminary data on the structure of the orthorhombic phase for YCUO2.5 suggests 
two sets of locally triangular configuration of copper sites with 2 ~ ions in the center 
of one of the sets of triangles, providing superexchange paths (Fig.l). The two sets of 
triangles are separated and have no 02— ions in their midst, so we have a good possibility 
of one-dimensional exchange coupled Cu 2+ i.e. a S = 1/2 system. We have carried out 
NMR measurements of the relaxation rates 1/Ti, 1/T 2 between 70° K and 230°-ftT. It is 
evident from the NMR results that there are substantial exchange couplings between the 
S = 1/2 Cu 2+ spin moments in this structure. We find an activated behavior as a function 
of T with a gap ~ 650° K. Such couplings would arise between nearest neighbor spins to an 
2 ~ ion, which would act as the conduit for 120° exchange paths. If we presume that such 
couplings between nearest neighbor and second-neighbor spins to an 2 ~ ion are negligible 
(60° exchange paths), then the system divides into a series of parallel one-dimensional 
"sawtooth" lattices of exchange coupled spins, which are only weakly interacting. Such a 
system, the sawtooth lattice is analyzed in this paper, and shown to have an interesting 
and unique set of magnetic excitations, namely quantum solitons very similar to the ones 
discussed in SS, with a new feature, namely the kink antikink symmetry in the Majumdar 
model is broken here. 



II. THE SAWTOOTH CHAIN: ANALYTICAL RESULTS 



The Hamiltonian for the sawtooth chain may be written as a sum of Hamiltonians 
governing triangles of spins (see Fig. 2) 

N 

H = H n 

n=l (2.1) 
H n = — ^ ( §2n-l + S 2n + S 2n +1 ) 2 — - ] ■ 



Here TV is the number of triangles and J denotes the antiferromagnetic coupling. We may 
consider either open chains (with 2N + 1 sites) or periodic chains (with 2N sites). 

To find the ground states of (2.1), we note that H n is proportional to a projection 
operator since its eigenvalues are (if the total spin in triangle n is 1/2) and 3 J/2 (if 
the total spin is 3/2). For two sites i and j with i < j, we denote the singlet state 
mathematically by = {\aij3j) — |/3jOj))/-\/2, and pictorially by a double line joining i 
and j as indicated in Figs. 3 and 4. (Here and (3i denote spin up and down respectively 
at site i). Then the states with total spin 1/2 in triangle n can be thought of as either 
[2n — 1, 2n] with the spin at site 2n + 1 free (either a or or as [2n, 2n + 1] with the spin 
at 2n — 1 free. The other possible pairing [2n — 1, 2n + 1] is linearly dependent since 

[2n, 2n + 1] a 2n = [2n - 1, 2n] a 2n +i + [2ra, 2n + 1] a 2n -i (2.2) 

It is now easy to show that the periodic chain has two degenerate ground states (with 
energy E Q = 0) given by 

JV 

| I ) = Y[ [2n-l,2n] 

n=1 (2.3) 
and | i7 ) = JJ [2n, 2n + 1] , 

n=l 

where SW+i = [4, 8]. The open chain has 2(N + l)degenerate ground states (with 
E a = 0) given by 

m N 

| 2n + 1, a or /3 ) = ( ][| [2n - 1, 2n] ) ( ]~[ [2n, 2n + 1] ) times 

n=l n=m+l (2-4) 

a2m+i or /3 2 m+1 ? 

where m may take any value from to N. Such a state is shown in Fig. 3. This 
configuration can be thought of as a kink at site 2m + 1 which separates the ground state 
I on its left from the ground state II on its right. Since all states in (2.4) have the same 
energy (namely, zero), a linear combination of them like 

1 N 

| k ) = —j= exp(z/cm) | 2m + 1, a ) (2-5) 

V m=0 



has the same energy for all values of k. In the limit N — > oo, (2.5) denotes a momentum 
eigenstate. We therefore see that kinks in the sawtooth chain have the dispersionless 
spectrum u(k) = 0. Further, kinks only differ from ground states I and II at a single site. 

We will see below that antikinks are quite different in that they have a non-trivial 
dispersion and they do not just differ from the states I and II at only one site. In fact, 
we are unable to solve for the wave function and dispersion of antikinks exactly. The 
best we can do is to variationally estimate these quantities more and more accurately by 
considering antikink configurations spread over 1 site, 5 sites and so on. 

We would like to make a few comments before examining the antikinks. Firstly, it can 
be shown by induction that the states in (2.3 - 2.4) are indeed the only ground states for the 
sawtooth chain [9] . Secondly, it can be proved rigorously that there is a finite gap between 
the ground states and the first excited state [10, 11]. Our discussion of antikinks will lead 
to an accurate estimate of this gap. Finally, the situation here may be contrasted with 
that obtaining in the railroad trestle which was first studied by Majumdar and coworkers 
[5]. The Hamiltonian for the railroad trestle differs from (2.1) in also having a H n for the 
sites (2n, 2n + 1, 2n + 2). As a result, this model only has the two ground states of types 
I and II (except possibly for free spins at the ends if the chain is open). There is a finite 
gap to excited states. Kinks and antikinks are on the same footing in the railroad trestle. 
They are not exactly solvable but they can be shown to have identical dispersions. 

We now study antikinks to a first approximaton by considering a state like the one 
shown in Fig. 4 (a). This is a configuration which has ground state II to the left of a 
1-site cluster (located at site 2n) and ground state I to its right. We denote this state by 
| 2n ) i . (For simplicity of notation, we will henceforth drop the spin index, a or /?, of the 
free spin). In the limit N — > oo, we consider a momentum eigenstate 



k) 



1 



n 



exp(z/cn) \ 2n)i . 



(2.6) 



This state has a non-trivial norm because 



!< 2n \ 2m) l = {-l) n ~ m / 2^~ m \ . 



(2.7) 
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Hence 

(k\k) = 3/(5 + 4 cos k) . (2.8) 

Further 

1 (2n\H l \2m) 1 = ^ J 8 nl 8 lm (2.9) 

implies that 

(k\H \k) = ^ J . (2.10) 
Our estimate of the dispersion based on this 1-cluster calculation is therefore 

u(k) = (- + cos/c) J . (2.11) 

(The kinks and antikinks in the railroad trestle have the same dispersion as in (2.11) for 
the 1-cluster approximation). Eq. (2.11) has a minimum at k = it where the antikink 
energy is J/4. This is our first estimate of the gap in the open chain and, as we will argue 
below, in the periodic chain also. 

We may now improve our estimate by considering 3-cluster configurations. Due to 
Eq. (2.2), however, a cluster of three neighboring sites with 5 = 1/2 can be reduced to 
a superposition of 1-cluster states like | 2n )i . So we have to continue on to 5-cluster 
configurations. The only two linearly independent configurations that we need to consider 
are the ones shown in Figs. 4 (b) and (c). We denote these two by | 2n ) 2 and | 2n ) 3 
respectively where 2n denotes the center of the 5-clusters. Note that | 2n ) 2 and | 2n ) 3 
are related to each other by reflection about the site 2n. We now consider a momentum 
eigenstate with two complex variational parameters a and b 

| k,a,b) = — L V" exp(ikn) | 2n ) x + a | 2n ) 2 + b | 2n ) 3 , (2.12) 

and minimize u(k,a,b) by varying a and b. The computation is straightforward though 
lengthy. We first obtain the overlaps 

2 (2n\2m) 1 = / 2l n " m l if n>m + l, 

= -1/4 if n = m , (2.13) 
= - (-l) n " m / 2l n " m l if n<m-l, 
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2 (2n | 2m ) 2 = 1 if n = m , 



= —1/8 if |n — m| = 1 , 

= -\ (-1)™ I 2l n " m l if |n-m|>2, 



(2.14) 



2 (2n | 2m ) 3 = (-l) n " m / 2l n " m| if n > m + 2 , 

= —1/8 if n = m+l, 

= -1/2 if n = m , (2.15) 

= 1/4 if n = m — 1 , 

= ^ (-l) n " m / 2l n " m l if n<m-2. 



All other overlaps can be obtained by using the reflection symmetry mentioned above. 
Thus, 



3 ( 2n | 2m )i = 2 ( 2m | 2n )i, 
and s( 2n \ 2m ) 3 = 2 ( 2n | 2m ) 2 , 



(2.16) 



We therefore obtain 



(k\k) 
where Ai 



1 



(5 + 4cosA;)A 1 + (4 + 2e~ lk ) A 2 + (4 + 2e lk )A 



Ak 



5 + 4 cos k - 

1 - ^(a + a* + b + b*) + ^(aa* + bb*) (4-cosA;) 

- hab* + ba*) + lab*(2e- lk - e* fc ) + ^6a*(2e* fc - e" lfc 
2 8 8 



and A 2 = e — (-2 - a* + 2a - 6 + 26*) + %— 
4 16 



-2aa* - 266* + 4a6* + 6a*) . 

(2.17) 



Next we calculate the matrix elements of H n . Thus 



2 {2n- 2 | H n 


2n )i = 


8 


2 (2n + 2 | H n 


2n )i = 


3 -j 

8 


2 (2n-4 | H n -! 


2n ) 2 = 


16 


o(2n 1 1 


2n )? = 


- j 

4 


3 (2n-4 | tf n _! 


2n ) 2 = 




3 (2n | F n _i 


2n ) 2 = 


"I"' 


2 (2n-4 | # n _! 


2n> 3 = 





(2.18) 



All other matrix elements can either be obtained from the above by translation or reflection 
symmetry, or are zero. We then find that 



3 3 

( k I H I k ) = - + i - sink (a- a* - b + b* ) 
* 1 1 ' 4 4 

3 3 3 3 

+ (- - - cos 2k) (aa* + bb* ) + (- cos 2k - -)(ab* + ba* ) . 
2 8 8 4 



(2.19) 



We have found that the minimum value of u(k) occurs at k = it if we take a = b to 
be real. Then 

1 1 4- On 2 

(2.20) 



, N 1 l + 2o 2 
w(7r;a) = - 



4 1 -a + a 2 /2 ' 

This has a minimum at a = (3 — Vl~7) /4 = —0.2808 where u = 0.2192 J. (For the railroad 
trestle, (k\k) is the same as in (2.17) while (k\H\k) has the extra term 3(aa* +bb* )/4 on 
the right hand side of (2.19). Hence 

1 1 + 4a 2 



lo(tt 7 a) 



(2.21) 



4 1 -a + a 2 /2 ' 

whose minimum value is 0.2344 J ). 

We see that the estimate of the gap changes relatively little on going from 1-clusters to 
5-clusters. This is because of the small correlation length £ in this system. We expect that 
the estimate of the gap A(l) from an /-cluster calculation will differ from the true gap A(oo) 
by terms of order exp(— Z/£). The gap in the railroad trestle chain has been estimated from 



a 9-cluster calculation in Ref. 12. From the values A(l) = 0.25J, A(5) = 0.2344J and 
A (9) = 0.2340 J, we estimate that £ ~ 1.1. While we have not computed A (9) for the 
sawtooth chain, we expect that it will differ very little from A(5) for a similar reason. 

To summarize so far, we have seen that kinks have zero dispersion while antikinks have 
a dispersion with a gap of 0.2192 J at k = n. We now identify the latter figure with the 
gap in the open chain. This assumes that there are no bound states of several kinks and 
antikinks which have a lower energy. For the railroad trestle, it is in fact known that there 
is no bound state of a kink and an antikink which has lower energy than a well-separated 
kink and antikink [4]. 

We will now argue that a periodic chain has the same gap and, further, that it has 
a dispersionless spectrum for its lowest excitation. Any excitation in a periodic chain 
must consist of a succession of alternating kinks and antikinks. In the absence of low- 
energy bound states, the lowest excitation in a long periodic chain will consist of one kink 
well-separated from one antikink. The energy of this state is the sum 

u(Q) = Wif(fei) + u^(k 2 ) , (2.22) 

where k\ and k 2 denote the momenta of the kink (K) and the antikink (K), and the total 
momentum of this state is Q = ki + k 2 ■ It is now clear that since uixiki) = for all fci, 
the minimum possible value of lo(Q) is given by the constant A = uj-^ (n) since we can 
always choose ki = Q — n. 

Indeed, numerical studies of finite periodic chains upto N = 10 by Kubo had indicated 
the existence of a dispersionless spectrum with uj(Q) ~ 0.219 J for all Q [13]. We now 
have the explanation of this striking property of the periodic sawtooth as arising from 
the dispersionless spectrum of the kink. Further, our 5-cluster computation has already 
yielded an estimate of the gap which is very close to the value obtained numerically. 

We may now use the above results to study low-temperature thermodynamic proper- 
ties of the sawtooth chain. For instance, we can estimate the magnetic susceptibility based 
on the picture of a low density of alternating kinks and antikinks which are well-separated 
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and noninteracting. In the presence of an external magnetic field B, they have the ener- 
gies —2[iS z and toj^(k) — 2fxS z respectively where \x is the Bohr magneton. If we use the 
1-cluster expression for Wj^(k) given in Eq. (2.11), we obtain the partition functions for 
one kink and one antikink as 

xk = 2 cosh( / u/3S) 

5 (2-23) 
and x-^ = 2 cosh(/i/3£?) I {(3J) exp( — -^(3 J) 

respectively, where (3 = 1/ksT is the inverse temperature and I Q is a modified Bessel 
function. The free energy per site is then given by 

/ = - j (x K x« ) 1 / 2 , (2.24) 

and the magnetic susceptibility per site x = —(9 2 f / dB 2 )b=o is 



xJ_ 

2/x 2 
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1/2 



J (/3J)exp(--/3J) . (2.25) 



Note that this thermodynamic quantity exhibits a gap equal to J/8 at very low temperature 
which is half the sum of the gaps for the kink (zero) and the antikink (J/4). In Fig. 5, we 
show the magnetic susceptibility as a function of ksT/ J. 

Coming back to the system YCu0 2 . 5 , we see that the gap of ~ 0.22J, if equated to 
the NMR activation energy, implies that J ~ 3000^, which is rather too large. Indeed, 
the largest J's are a factor of 2 smaller than this, as measured in the high T c systems, 
which have comparable Cu — O bond lengths as in these compounds, namely ~ 2A° . It 
then seems likely that these systems either do not allow for a decoupling between these 
sawtooth lattices, forcing say a pair of these excitations, or else there might be pairwise 
dimerization, which could be signalled in detailed structural studies. It would be interesting 
to study these issues in detail experimentally, as would the ESR on these systems at 
elevated temperatures, say 700°K, where the four-fold degeneracy would be reflected in 
the existence of free spin 1/2 excitations with a thermally activated density. 
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FIGURE CAPTIONS 



1. A picture of the CuO planes in YCUO2.5 from preliminary structural data. 

2. The sawtooth chain. The three sites forming triangle n are numbered as shown. 

3. A kink configuration with the free spin at site 2n + 1. The double lines indicate 
singlets. 

4. The three antikink configurations centred about site 2n. (a) is a 1-cluster while (b) 
and (c) are 5-clusters. 

5. Low-temperature magnetic susceptibility as a function of UbT/J. 
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Fig. 2 
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